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ABSTRACT 

Gamma-ray burst sources are distributed with a high level of isotropy, which 
is compatible with either a cosmological origin or an extended galactic halo 
origin. The brightness distribution is another indicator used to characterize 
the spatial distribution in distance. In this paper we discuss detailed fits of 
the BATSE gamma-ray burst peak-flux distributions with Friedmann models 
taking into account possible density evolution and standard candle luminosity 
functions. A chi-square analysis is used to estimate the goodness of the fits and 
we derive the significance level of limits on the density evolution and luminosity 
function parameters. Cosmological models provide a good fit over a range of 
parameter space which is physically reasonable. 


1. INTRODUCTION 

The gamma-ray bursts are distributed isotropically on the sky (Fishman et 
al. 1994; Tegmark et al. 1996; Briggs et al. 1996; although a recent paper shows 
that at least one quadrupole term is non zero, Balazs, Meszaros and Horvath 
1998), which is compatible with either a cosmological origin or an extended 
galactic halo origin; where the first origin seem to be more probable (Tegmark 
et al. 1996, Briggs et al. 1996). Recently, the successful identifications made 
by the Beppo-SAX satellite, followed by the detection of optical counterparts 
(van Paradijs et al. 1997), seem to give a firm support for the models aiming to 
explain the bursts by merging neutron stars (Usov and Chibisov 1975; Rees and 
Meszaros 1994; Meszaros and Rees 1997) and seem to put them definitely into 
the cosmological distances. In addition, the fainter bursts seem to have longer 
durations (Norris et al. 1995; Meszaros et al. 1996). Hence, all these seem to 
put the bursts definitely into the cosmological distances. Nevertheless, further 
proofs of cosmological origin are still useful. The search for a further support is 
the aim of this paper. 

The brightness distribution is the indicator used to characterize the spatial 
distribution in distance, and this can also be used to test the cosmological 
hypothesis. This is generally done by investigation of the functional behavior of 
the integral number N of sources with peak photon flux rates P above a certain 
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value, N(> P). Comparisons of observed versus expected values in Friedman 
cosmologies have been discussed, e.g., by Mao and Paczynski (1992), Dermer 
(1992), Piran (1992) and Wasserman (1992). Statistical fits to a logiV — logP 
distribution have been done by Wickramasinghe, et al. (1993), Cohen and 
Piran (1994), Ernslie and Horack (1994), Meszaros and Meszaros (1995), Loredo 
and Wasserman (1995), Horack, Ernslie and Hartmann (1995), Fenimore and 
Bloom (1995), Horvath, Meszaros and Meszaros (1996), Meszaros and Meszaros 
(1996). One of the main questions that such fits must address is the size of the 
parameter space region which is compatible with a cosmological distribution, 
and whether such parameters are reasonable. If the acceptable region contains 
physically plausible parameters and is not too restricted, one may assume the 
consistency of the observations with a general type of models; if, on the other 
hand, the acceptable region is very small and/or populated mainly by physically 
implausible parameters, fine-tuning would be required to fit the observations, 
and the case for consistency with those models is weaker. Such consistency, and 
absence of fine-tuning, is a requirement expected of any successful model of the 
GRB distribution, whether cosmological or galactic. Here, we shall address only 
the question of the consistency of the number distribution under the hypothesis 
of a cosmological distribution with a standard candle luminosity distribution. 

Most cosmological fits have been made with relatively specialized models, 
generally either with non-evolving or evolving density standard candle models, 
or with non-evolving luminosity functions. Limits on the luminosity function 
were investigated in cosmology with a pure density evolution by Horack, Ernslie 
and Hartmann (1995) using a method of moments. In Euclidean space, lim¬ 
its have been investigated by Horack, Ernslie and Meegan (1994), Ulmers and 
Wijers (1995) and Ulmer, Wijers and Fenimore (1995). Most cosmological cal¬ 
culations have used either the IB or the 2B BATSE data base. In the present 
paper we make detailed chi-squared fits of the observed brightness distribution. 
This is done both using the BATSE 2B catalogue of sources (Meegan et al. 
1996), and combining the BATSE catalogue with information published for the 
PVO counts (Fenimore and Bloom 1995). The significance levels of the various 
cosmological fits is discussed for both the 2B and the expanded burst sample. 

Note that, as it seems, there are several subclasses of gamma-ray bursts 
(Kouveliotou et al. 1993; Dezalay et al. 1996; Pendleton et al. 1997; Meszaros 
et al. 1998; Bagoly et al. 1998). Nevertheless, the criteria are still ambiguous. 
Therefore, in this article these subclasses will not be considered. 


2. MODELS AND FITS 

Analytical expressions for the integral burst number counts N(> P) with 
peak photon flux rate in excess of P (units of photon cm -2 s -1 ) were discussed 
by Meszaros and Meszaros (1995) and by Meszaros and Meszaros (1996) for 
arbitrary Friedmann models with zero cosmological constant . As discussed in 
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Meszaros and Meszaros (1995), effects of a non-flat cosmology (0 o < 1) are 
small, and to a first approximation can be neglected. Below we assume Q 0 = 1 
everywhere. The effect of a pure density evolution is approximated through a 
dependence 

n(z) = n 0 (l + z) D , (1) 

where n is the physical burst density rate in cm^ a yr~ 1 , n 0 is the density rate 
at z = 0 ( D = 3 corresponds therefore to a non-evolving, constant comoving 
density). It is assumed that the sources have the same intrinsic luminosity C Q . 
Therefore the photon luminosity function in the 50-300 keV range is represented 
by the form, 

$(£) = n 0 S(C - Co) (2) 

For the rest of the paper we take P to be peak photon flux [c?n -2 s -1 ], n 0 is the 
physical density of bursts per year at 2 = 0. 

The integral number distribution of bursts per year with peak flux rate above 
P is given by Meszaros and Meszaros (1995) as 


N(> P ) 


47t Cl /2 n 0 
3 (4ttP) 3 / 2 


( 3 ) 


where / is a dimensionless analytical function of P and the model parameters, 
i.e. the luminosity function parameters C 0l the density evolution parameter D, 
and the density ?r(0) = n 0 . 

For the numerical fits we used the BATSE 2B catalog (Meegan et al. 1996). 
The 1024 ms peak fluxes P (photons cm~ 2 s~ 1 ) were used, and only events 
with peak count rates divided by threshold count rates C max /C m i n > 1 were 
included, where C m i n is the published count threshold for each event (Mee¬ 
gan et al. 1996). The 2B sample with this criterion consists of 278 entries in 
the catalog. Applying the efficiency tables published with the catalog to cor¬ 
rect for detector inefficiency near the trigger threshold, the nominal number 
of bursts accumulated by BATSE over a period of two years with peak fluxes 
above logP > —0.6 is 369. We chose for these bursts a binning equidistant in 
log 10 P, with step size 0.2 between -0.6 and 1.2, which gives 9 equal bins with 
a minimum number of 7 events per bin (in the highest P bin, log 10 P = 1.0 to 
1.2) for the two year 2B sample. The fits were made to the differential burst 
number distribution V(P) as a function of peak photon flux P (since only in the 
differential distribution may the bins be considered independent of each other 
for a x 2 fit) and the errors in each bin were taken to be the square root of the 
number of events in that bin. 

Some of the fits were made using an extended 2B plus PVO sample. For 
the PVO events, we used the PVO portion of Table 2 of Fenimore and Bloom 
(1995) for log 10 P > 1.2. A number of subtle issues concerning a matching 
between the different PVO and BATSE data sets are discussed by Fenimore et 
al. (1993), who indicate that systematic uncertainties of ±10% in the relative 
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normalization cannot be ruled out. The matching of the level of the BATSE 
and PVO curves was taken directly from Fenimore and Bloom (1995). The 
PVO data was rebinned, ignoring PVO bursts below log 10 P = 1.2 so as not 
to count twice, and its level was renormalized so that the matching 2B data 
had the same level as in the original 2B catalog, i.e about 2 years. The errors 
for the PVO sample were also renormalized taking into account the fact that 
data had accumulated over more than ten years in the PVO case, keeping the 
relative errors the same. We used 5 bins in the PVO range, so that the combined 
2B+PVO fits have 9+5=14 bins, reaching up to log 10 P = 3.0. 

The fits (standard candle with density evolution) involve the fewest parame¬ 
ters: the photon luminosity C 0 (ph/s), the density n 0 and the density evolution 
parameter D. For the 2B sample between peak fluxes —0.6 < logP < 1.2, 
the free parameters are p = 3 the degrees of freedom are / = 6 and the best 
Xred, (reduced chi-square or y 2 divided by degrees of freedom) is 0.85 at the 
innermost mark. The lcr, 2cr, 3cr significance contours were determined using 
the standard prescription (e.g., Press, et ah, 1986, or Lampton, Margon and 
Bowyer 1976). The fit (Fig. a) is good over an elongated region describing a 
relation between the luminosity and the density evolution. For faster density 
evolution n oc (1 + z) D (larger D) the luminosity must increase because most 
sources are at larger redshifts, while for slower or negative evolution the lumi¬ 
nosity must decrease, since most sources are at small redsliift {D = 3 is constant 
comoving density). The optimal fit is obtained for D = 3.5 and C 0 ~ 10 57 s -1 . 
This luminosity is close to the standard candle value deduced, e.g. Horack, 
Emslie and Hartmann (1995) and Fenimore and Bloom (1995). (corresponding 
to C 0 ~ 10 51 er( 7 s -1 for typical photon energies of 0.5 MeV). However the pref¬ 
erence for D — 3.5 was not, as far as we can tell, encountered in previous fits. 
The Xred around the best fit minimum is 0.85; however, the lcr region around 
it is rather large, even if not very wide, so this preference is not strong. 

The fits using the 2B+PVO sample (p = 3, / = 11) are shown in Fig. b, with 
a best xled = 0-62 at the central mark enclosed by its lcr, 2a, 3cr contours. We 
note that this ignores any possible systematic errors in matching BATSE and 
PVO beyond what is done in Fenimore et al. (1993), and Fenimore and Bloom 
(1995). If any extra errors were present, they could in principle increase the 
size of the significance regions discussed below (e.g. it might add an extra free 
parameter for the relative normalization). However, such errors are extremely 
difficult to quantify without going into additional details of the instruments, and 
we follow Fenimore and Bloom (1995) in adopting their relative normalization 
as adequate without further manipulation. The effect of the rare high flux 
PVO bursts satisfying a tight N oc P~ 3 / 2 correlation at 1.2 < log 10 P < 3.0 
is to improve the best fit (lower Xred) an d to place it at a somewhat smaller 
luminosity £„~ 5 x 10 56 s^ 1 and closer to comoving constant density evolution, 
D ~ 3. This is in good agreement with Fenimore and Bloom’s (1995) value 
of L 0 ~ 5 x 10 50 ergs -1 . However, as seen from Fig. b, the lcr region around 
this best fit minimum is compatible with both larger and smaller C 0 and D. In 
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contrast to the pure 2B fit, however, the joint lcr upper limit for C 0 and D are 
C 0 < 5 x 10 57 s -1 , D < 4.5 (or 3 ct joint upper limits C Q < 5 x 10 58 s -1 , D < 5). 


3. CONCLUSION 

The fits presented above show that a cosmological interpretation is compat¬ 
ible with the data under a variety of assumptions. Good fits of the observed 
differential distribution of bursts N(P ) as a function of peak photon flux P are 
obtained under a standard candle luminosity function assumption. Fits were 
obtained for a range of density evolution indices D 1 defined through a physical 
density dependence n 0 oc (1 + z) D where D = 3 is equivalent to a non-evolving, 
constant comoving density case. The cosmological fits obtained are of good 
quality (Xred ~ 1) f° r a ran ge of plausible model assumptions. 

The present results are in significant agreement with the literature (Horack, 
Emslie and Meegan 1994, Horack, Emslie and Hartmann 1995, Ulmer, Wijers 
and Fenimore 1995, Hakkila et al. 1995, Horvath, Meszaros and Meszaros 1996). 
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Figure Caption 

Cosmological fits for f2 0 = 1, A = 0 for standard candle photon luminosity 
C 0 (photon/s) and physical density evolution n(z) oc (1 +z) D . The inner dots are 
the best fit Xred minimum location, with la, 2 a, 3<r contours increasing outwards, 
a) Top: using the BATSE 2B data base; b) Bottom: using the 2B plus PVO 
information. 
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